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Abstract 

We analyze the control of the propagation speed of three wave packets interacting in a medium 
with quadratic nonlinearity and dispersion. We found analytical expressions for mutually trapped 
pulses with a common velocity in the form of a three-parameter family of solutions of the three- 
wave resonant interaction. The stability of these novel parametric solitons is simply related to the 
value of their common group velocity. 

PACS numbers: 05.45.Yv, 42.65.-k, 42.65. Sf, 42.65. Tg, 52.35.Mw 
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A three-wave resonant interaction (TWRI) is ubiquitous in various branches of science, as 
it describes the mixing of waves with different frequencies in weakly nonlinear and dispersive 
media. Indeed, TWRI occurs whenever the nonlinear waves can be considered as a first- 
order perturbation to the linear solutions of the propagation equation. TWRI has been 
extensively studied alongside with the development of nonlinear optics, since it applies to 
parametric amplification, frequency conversion, stimulated Raman and Brillouin scattering. 
In the context of plasma physics, TWRI describes laser-plasma interactions, radio frequency 
heating, and plasma instabilities. Other important domains of application of TWRI are 
light-acoustic interactions, interactions of water waves, and wave-wave scattering in solid 
state physics. Two classes of analytical soliton solutions of the TWRI have been known for 



over three decades. The first type of solitons describes the mixing of three pulses w 
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with their respective linear group velocity, and interact for just a short time 
second type of solitons, also known as simultons, are formed as a result of the mutual trapping 
of pulse envelopes at the three different frequencies. Hence the three wave packets travel 
locked together with a common group velocity j^. In all of the above discussed domains 
of application, parametric TWRI solitons play a pivotal role because of their particle-like 
behaviour, which enables the coherent energy transport and processing of short wave packets 

].QQ 

In this Letter we reveal that the class of TWRI simultons (TWRIS) is far wider than 
previously known. We found a whole new family of bright-bright-dark triplets that travel 
with a common, locked velocity. The most remarkable physical property of the present 
solitons is that their speed can be continuously varied by means of adjusting the energy of 
the two bright pulses. We studied the propagation stability of TWRIS and found that a 
stable triplet loses its stability as soon as its velocity decreases below a well defined critical 
value. Another striking feature of a TWRIS is that an unstable triplet decays into a stable 
one through the emission of a pulse, followed by acceleration up to reaching a stable velocity. 

The coup,ed partia, differential equations (PDEs) representing TWRI in (1 + 1) dimen- 
sions read as |2|: 

En — VxE lz = E^E^, 

E-n — ViE^z = —E{El, (1) 
E3 t — V 3 E 3z = E^E^, 



2 



where the subscripts t and z denote derivatives in the longitudinal and transverse dimension, 
respectively. Moreover, E n = E n (z,t) are the complex amplitudes of the three waves, V n are 
their linear velocities, and n = 1,2, 3. We assume here the ordering V\ > V 2 > V3 which, 
together with the above choice of signs before the quadratic terms, entails the non-explosive 
character of the interaction. In the following, with no loss of generality, we shall write the 
equations Q in a reference frame such that V 3 = 0. A remarkable property of the Eqs.Q 
is their invariance with respect to the transformation 

E n (z, t) = s exp[i(q n z n + a n )]E n (sz + z , st + t ) (2) 

where a\+a 2 + a 3 = 0, q n = q(V n +i — V n+ 2), z n = z + V n t are the characteristic coordinates 
and n = 1, 2, 3 mod(3). As the transformation (J2J depends on six real parameters, namely 
ai,a 2 , s,q, z and t , clearly one may introduce these parameters in the expression of any 
given solution of the TWRI equation. 

The evolution equations represent an infinite-dimensional Hamiltonian dynamical 
system, with the conserved Hamiltonian 

1 



H 

4i 



dz, 



+^3(E^ Z E 3 — E 3z E^) + 2EiE 2 E 3 — 2E\E* 2 E. 
energies (Manley-Rowe invariants) 



(3) 
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1x2 = h + h = 5 / (i^il + \Ei\ l )dz, (4) 



00 



1 r+00 

l2 3 = h + h = - (\E 2 \ 2 + \E 3 \ 2 )dz, (5) 



and total transverse momentum 



^ r+00 

J = — J [{E{Eu - EiE* lz ) - (E 2 E 2z - E 2 E 2z 
+ (E* 3 E 3z - E 3 E* 3z )}dz. 



(6) 



Each of the above conserved quantities is related to a given internal parameter of the TWRIS 
which, in turn, is associated with a symmetry (e.g., phase rotation or space translation) of 
the TWRI equations ((TJ) 7]. As a consequence, one may expect that Eqs.JTJ) possess a 
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three-parameter family of soliton solutions. We found such soliton solutions by using the 
recent results on the TWRI equations as presented in Ref.0]. Their expression is 

1 y/W+W^Vi - v 2 ) cosh[B{z + vt)} ' 1 a) 

2 y/\b\ 2 + \a\ 2 g{Vi-V 2 ) cosh[B{z + Vt)} ' 1 ' 
= {1 + T^^ 1 - ^h[B(z + Vt)}}} ^ _ - } (7c) 

where 

b=(Q- l){p + ik/Q), r = p 2 -k 2 -\a\ 2 , 



Q = l\ll [r + v^+^V]> 
E = + y 2 _ Q(y : _ v 2 ) }/(Vx - V 2 ), 

V = 2V 1 V 2 /[V 1 + V 2 - Q(Vx - V 2 ) }, 
X = k[V 1 + V 2 - (V, - V 2 )/Q\/{V X - V 2 ), 
lu = -2W x V 2 jiy x - V 2 ), q n = q(V n+1 - V n+2 ), 
9n = \iy n - V n+1 ) (V n - V n+2 )\' 1/2 , g = gi g 2 g 3 . 



and n = 1,2,3 mod (3). The above expressions depend on the five real parameters 
Vi,V 2 ,p,k,q, and the complex parameter a. From the definition of Q, one can see that 
the above parameters must be chosen so that if k = 0, then p 2 > \a\ 2 . 

The TWRIS is composed of two bright pulses (|7al I7b|) . and a kink or shock-like pulse 
(|7cj). which travel with a common locked velocity V. The expressions flJJ) may be represented 
in a more convenient form as 

£n(£,r) = £/„(£) exp[z$ n (£,T)], n = 1,2,3. (8) 

Here we used a reference frame that moves along with the soliton, with co-ordinates £ = 
z + Vt, r = t where U and <3> are real functions and $ n (£, t) = 4> n r +/ n (£). A simple analysis 
of ((HI) shows that, for any value of the parameters, the pulse amplitudes U x (£), ^(0 an d 
C/ 3 (£) are even functions of £ and the phase constants satisfy 4>i + 4> 2 + (fi 3 = 0. On the 
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FIG. 1: (color online). Amplitude and phase of two simultons (|7|) at t = with V\ = 2, V% = 1, 
a = 1, k = 0.5, q = 0.5; in the upper figures p = 1, in the lower figures p = —1. Thin solid curve 
E\, dashed curve E 2 , thick solid curve E3. 

other hand, if k — the phase profiles are all piecewise linear in £ and obey the condition 
/i(£) + MO + /a(0 = for p£ > and A(0 + / 2 (£) + / 3 (£) = tt for p£ < 0. Whereas for 
k 7^ the phase profile /3(C) is nonlinear and cos[/i(£) + f 2 (0 + /3(C)] * s an °dd function of 
£; moreover the kink pulse S3 is "grey" if k 7^ and is "dark" if k = 0. Such amplitude and 
phase front profiles prevent a net energy exchange among the three waves. It is important 
to point out that the condition — 1 < Q < 1 leads to a speed V that lies in-between the 
characteristic velocities V\ and V 2 of the two bright pulses, i.e. V\ > V > V 2 . The above 
described properties mean that TWRIS represent a significant generalization with respect 



to previously known three-wave solitons which ex 



libit a simple (constant) phase profile and 
j|. In FigHJwe plotted two characteristic 



correspond to the special case fc = g = 0,r>0 
examples of TWRIS amplitude and phase-fronts (0) . 

It is interesting to consider the physical meaning of the various TWRIS parameters ap- 
pearing in (JJJ). For a given choice of the characteristic linear velocities V\ and V2, we 
are left with the four independent parameters p, k, q, and a. We may note that p is ba- 
sically associated with the scaling of the wave amplitudes, as well as of the coordinates 
z and t. The parameter a determines the amplitude of the asymptotic plateau of the 
kink E3. The value of k provides the wave-number of a "carrier-wave". The parameter 
q adds a phase contribution which is linear in z and t. Since the system (JIJ) is invariant 
under a transformation (j2J), without loss of generality we may set a = 1, which reduces 
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FIG. 2: (color online), (a) phase constant 0i versus velocity V for different energies 1\ 2 . (b) energies 
1%, I 2 versus velocity V for representative values of the phase constant 4>\. In all cases <p 2 = —<fii, 
V x = 2,V 2 = 1. 

the number of essential parameters to just three, corresponding to the three symmetries 
of Eqs.Q. The parameters p,k,q in (JJJ) may be more conveniently mapped into the pa- 
rameters V,0i,02 of Eq. (JH|) . which provide a more direct physical insight into the features 
of a TWRIS. Such a mapping is obtained by comparing Eqs.(J7|) with (JSJ), and reads as: 
V = 2V 1 V 2 /[V 1 + V 2 -Q{V 1 -V 2 )}, <p 1 = qV 2 (V 1 -V)+ X V + L J , <\> 2 = qV x {V -V 2 ) - X V-u. 
The TWRIS is thus simply expressed in terms of its velocity V and the two phase constants 
0! and 02- Let us investigate what are the TWRIS properties for a fixed choice of the linear 
velocities V\ and V 2 , upon variations of its energy flows and transverse momentum. As an 
example, Fig. EJa) shows the dependence of the phase constant <pi on the locked velocity 
V, for the case where <p 2 = — 0i, with different values of the conserved energy 1\ 2 . More- 
over, Fig. EJb) illustrates the dependence of the energies I\ and I 2 (which happen to be 
time-independent for a TWRIS) on the locked velocity V, for different choices of the phase 
constant <pi = —02- As it can be seen, the intensity and phase profiles, as well as the energy 
distribution among the different wave packets, strongly depend upon the value of the locked 
velocity V. 

The next crucial issue is the propagation stability of TWRIS. A first insight into this 
problem may be provided by performing a linear stability analysis (LSA) as in Ref. [9|. Let 
us consider a perturbed TWRIS of the form 

r) = (x n (0 + P n (£, r))e^ T , n = 1, 2, 3 

where x n (£) = U n (^) exp[i/ n (£)] is the soliton profile, and we consider a weak perturbation 
\P n \ « \x n \. By inserting the above ansatz in Eqs.Q, and by retaining only linear terms 
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FIG. 3: (color online). Propagation of stable and unstable simultons. The common parameters 
are V\ = 2, V2 = 1, a = 1, k = 0.5, q = 1. For the stable simulton (upper figures) p = 1 
(V = 1.68 > V cr , (pi = 0.5, 02 = 2.2), whereas for the unstable simulton (lower figures) p = — 1 
(V = 1.1 < y cr , 0i = 1.43, 02 = -0.33). 

in P n , one obtains a linear system of PDEs. For the numerical analysis, these PDEs can be 
reduced to a system of ordinary differential equations -P(t) = MP(t), by approximating the 
spatial derivatives with finite differences, where P is the perturbation vector sampled on a 
finite grid. A necessary condition for the instability of a stationary solution x n (£) is that the 
matrix M has at least one eigenvalue with positive real part. Numerical computations over 
a wide parameter range show that eigenvalues of M exist with a positive real part whenever 
p < 0. On the other hand, for p > the largest real part of the eigenvalues is equal to zero, 
which means that the TWRIS are only neutrally stable. Note that the instability condition 
p < leads to the inequality V < V cr = 2ViV2/(Vi + V2). Extensive numerical integrations 
of Eqs.Q confirm that TWRIS with V < V cr (V > V cr ) are always unstable (stable). The 
propagation of either stable or unstable TWRIS is illustrated in Fig. which shows the 
general feature of unstable solitons with V < V cr . Namely, the simulton decays into a stable 
soliton with V > V cr , and it emits a pulse in the wave E 3 . It is quite remarkable that 
the dynamics of the decay from unstable into stable solitons may be exactly described by 
analytical solutions with variable velocity or boomerons A complete description of the 
parametric boomerons will be the subject of a more extended report. 

We performed further investigations of TWRIS stability by carrying out a multi-scale 
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asymptotic analysis (MAA) 1(J, The MAA aims to find the rate of growth A (with 



A << 1) of small perturbations, by linearizing Eqs.Q around the soliton solution. This 
procedure leads to a linear eigenvalue problem, whose solution can be expressed as an 
asymptotic series in A [12]. In this way, one obtains the following condition which defines 
the borderline between stable and unstable TWRIS 

<9(/ 12 ,7 23 ,7) 

where G is the Jacobian of the constants of motion Ii 2 , J23, J with respect to <f>i, 2 , V. Note 
that in Q I23 and J are obtained by re-normalizing the divergent integrals (|HJ) and (JHJ) 
according to the prescription 

1 



23- 2 



+00 

(|^| 2 + |£7 3 r-|^o| 2 )^, (10) 



J = J- 




where I-E30I = lirn^oo \E^\ is the asymptotic amplitude of the kink |L2(. Note that the 
availability of exact soliton solutions allows for the analytical calculation of the above inte- 
grals, hence of the condition ©, which is an extension of the well-known Vakhitov-Kolokov 
criterion. Thus provides a sufficient stability condition, which can only be applied under 
specific constraints such as <f>\ = 2 . Indeed, in this case we find that the condition G = 
leads again to the previously found marginal stability condition V = V cr = 2ViV 2 / (Vi + V2). 



A direct insight into the global stability properties of TWRIS for all possible values of 
their parameters can be obtained by means of a geometrical approach Q| . Indeed, TWRIS 
may be obtained as solutions of the variational problem 

8{H + - (0x + 2 )7 23 - VI) = 0, (12) 

where 5 is the Frechet derivative. In other words, TWRIS represent the extrema of the 
Hamiltonian (JHJ), for a fixed value of the energies and momentum (V,4>i,4>2 represent La- 
grange multipliers). Stable triplets are obtained whenever such extrema coincide with a 
global minimum of H. Clearly, if multiple solutions exist with the same ii 2 , 1 23, J, the 
stable solution is obtained on the lower branch of H. In this framework the condition 
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FIG. 4: (color online), (a) Hamiltonian H versus energy lyi and (b) I\2 versus velocity V. Char- 
acteristic velocities are V\ = 2, V2 = 1 and <p2 = —0.2. Dashed (solid) curves, unstable (stable) 
solitons. 

(JHJ) corresponds to solitons such that the normal vector to the three-dimensional surface 
H = H(Ii2, 123, J) lies in the space H = const.. The above geometrical considerations 
permit the visualization of the stability boundaries when considering a projection of the 
hyper-surface H = H(Ii 2 , 123, J) on the plane (Ii2,H). For example, Fig. HI displays the 
dependence of H upon I 12 for the case (pi = —02, where the criterion Q cannot be applied, 
and in the case <p\ = 02- Here it is evident that the two branches of the Hamiltonian merge 
exactly at V — V cr : at this point, the normal to the H curve is also orthogonal to the 
vertical axis. Interestingly enough, Fig. 0] shows that the borderline TWRIS corresponds 
to a minimum of the bright pulses energy J 12 with respect to V . To summarize, we have 
shown that different numerical and analytical methods concur in predicting that the TWRIS 
stability is determined by the condition V > V cr . 

Let us briefly discuss the experimental conditions for the observation of TWRIS in nonlin- 
ear optics. For instance, when considering a three-wave oeo interaction in a 5 cm-long bulk 
PPLN sample with 18/im periodicity, the field envelope carriers Xei = 1.5//m, Xe2 = 0.8/xm, 
X E3 = 1.8/xm, pulse durations of about lps, TWRIS can be observed with field intensities 
of a few MW/cm 2 . 

In conclusion, we have described a novel three-parameter family of (1+1) bright-bright- 
dark soliton waves as exact solutions of the three-wave resonant interaction equation. These 
TWRI solitons exhibit nonlinear phase-fronts curvatures, and exist for a given range of their 
locked velocity and energy flows. Their propagation stability has been investigated with the 
upshot that stable triplets occur whenever their velocity V is greater than a certain critical 
value V cr . On the other hand, unstable solitons dynamically reshape into stable solitons with 
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higher velocity. The remarkable properties of these parametric solitons may open the way 
to new possibilities for the control of coherent energy transport in various physical settings. 
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